Rabi oscillations in a qubit coupled to a quantum two-level system 
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We consider the problem of a qubit driven by a harmonically oscillating external field while it is 
coupled to a quantum two-level system. We perform a systematic numerical analysis of the problem 
by varying the relevant parameters. The numerical calculations agree with the predictions of a 
simple intuitive picture, namely one that takes into consideration the four-level energy spectrum, 
the simple principles of Rabi oscillations and the basic effects of decoherence. Furthermore, they 
reveal a number of other interesting phenomena. We provide explanations for the various features 
that we observe in the numerical calculations and discuss how they can be used in experiment. 
In particular, we suggest an experimental procedure to characterize an environment of two-level 
systems. 



I. INTRODUCTION 

There have been remarkable advances in the field of 
superconductor-based quantum information processing 
in recent years p). Coherent oscillations and basic gate 
operations have been observed in systems of single qubits 

and two interacting qubits HSI1SIEHEI3- ° ne ° f 

the most important operations that are used in manipu- 
lating qubits is the application of an oscillating external 
field on resonance with the qubit to drive Rabi oscilla- 
tions 0- [J 0- 0- El ■ A closely related problem with great 
promise of possible applications is that of a q ubit coupled 
to a quantum harmonic-oscillator mode IT^. IT3. Il^j. 

Qubits are always coupled to uncontrollable degrees of 
freedom that cause decoherence in its dynamics. One 
generally thinks of the environment as slowly reducing 
the coherence of the qubit, typically as a monotonically 
decreasing decay function. In some recent experiments, 
however, oscillations in the qubit have been observed that 
imply it is strongly coupled to quantum degrees of free- 
dom with long decoherence times ^J, ^5|. The effects 
of those degrees of freedom have been successfully de- 
scribed by modelling them as quantum two-level systems 
(TLSs) [H H Mill. Since, as mentioned above, 
Rabi oscillations are a simple and powerful method to 
manipulate the quantum state of a qubit, it is important 
to understand the behaviour of a qubit that is driven 
on or close to resonance in the presence of such a TLS. 
Furthermore, we shall show below that driving the qubit 
close to resonance can be used to extract more parame- 
ters about an environment of TLSs than has been done 
in experiment so far. The results of this study are also 
relevant to the problem of Rabi oscillations in a qubit 
that is interacting with other surrounding qubits. 

Some theoretical treatments and analysis of special 
cases of the problem at hand were given in Refs. (l5Lll9| . 
In this paper we perform a more systematic analysis in 
order to reach a more complete understanding of this phe- 
nomenon. We shall present a few simple physical princi- 
ples that can be used to understand several aspects of the 



behaviour of this system with different possible choices 
of the relevant parameters. Those principles arc (1) the 
four-level energy spectrum of the qubit+TLS system, (2) 
the basic properties of the Rabi-oscillation dynamics and 
(3) the basic effects of decoherence. We shall then per- 
form numerical calculations that will agree with that in- 
tuitive picture and also will reveal other results that are 
more difficult to definitively predict otherwise. Finally, 
we suggest an experimental procedure where the driven 
qubit dynamics can be used to characterize the environ- 
ment of TLSs. 

The paper is organized as follows: in Sec. II we in- 
troduce the model system and the Hamiltonian that de- 
scribes it. In Sec. Ill we present a few simple arguments 
that will be used as a foundation for our numerical anal- 
ysis of Sec. IV, which will confirm that intuitive picture 
and reveal other less intuitively predictable results (note 
that a reader who is sufficiently familiar with the subject 
matter can skip Sec. III). In Sec. V we discuss how our 
results can be used in experiment. We finally conclude 
our discussion in Sec. VI. 



II. MODEL SYSTEM 

The model system that we shall study in this paper 
is composed of a harmonically driven qubit, a quantum 
TLS and their weakly coupled environment |2l| . We as- 
sume that the qubit and the TLS interact with their own 
(uncorrelated) environments that would cause decoher- 
ence even in the absence of qubit-TLS coupling. The 
Hamiltonian of the system is given by: 



H{t) = JT q (t) + #TLS + Hi + H; 



Env; 



(1) 



where H q and -Htls are the qubit and TLS Hamiltoni- 
ans, respectively; H\ describes the coupling between the 
qubit and the TLS, and H^nv describes all the degrees 
of freedom in the environment and their coupling to the 
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qubit and TLS. The (time-dependent) qubit Hamiltonian 
is given by: 

H q (t) = - ^4 q) - |^ q) (2) 
+ Fcos(u>t) (sin(9foi q) + cos (9 f o-( q) ) , 

where A q and e q are the adjustable static control param- 
eters of the qubit, a^P are the Pauli spin matrices of the 
qubit, F and u) are the amplitude (in energy units) and 
frequency, respectively, of the driving field, and Of is an 
angle that describes the orientation of the external field 
relative to the qubit <x z axis. Although we have used a 
rather general form to describe the coupling between the 
qubit and the driving field, we shall see in Sec. Ill that 
one only needs a single, easily extractable parameter to 
characterize the amplitude of the driving field. We as- 
sume that the TLS is not coupled to the external driving 
field, and its Hamiltonian is given by: 



TLS observable can be recast into the above form using 
a simple basis transformation, keeping in mind that such 
a basis transformation also changes the values of # q , #tls 
and 9f. 

We assume that all the coupling terms in H^nv are 
small enough that its effect on the dynamics of the 
qubit+TLS system can be treated within the framework 
of the markovian Bloch-Redfield master equation ap- 
proach. We shall use a noise power spectrum that can de- 
scribe both dephasing and relaxation with independently 
adjustable rates, and shall present our results in term of 
those decoherence rates. For definiteness in the numer- 
ical calculations, we take the coupling of the qubit and 
the TLS to their respective environments to be described 
by the operators ai a \ where a refers to the qubit and 
the TLS. Note, however, that since we use the relaxation 
and dephasing rates to quantify decoherence, our results 
are independent of the choice of system-environment cou- 
pling operators. 



_ AtlS (TLS) _ £TLS (TLS) 
TLS 2 <J X 2 O z 



(3) 



where the definition of the parameters and operators is 
similar to those of the qubit, except that the TLS param- 
eters are uncontrollable. Note that our assumption that 
the TLS is not coupled to the driving field can be valid 
even in cases where the physical nature of the TLS and 
the driving field leads to such coupling, since we gener- 
ally consider a microscopic TLS, rendering any coupling 
to the external field negligible. 

The energy splitting between the two quantum states 
of each subsystem, in the absence of coupling between 
them, is given by: 



E n 



(4) 



where the index a refers to either the qubit or the TLS. 
The corresponding ground and excited states are, respec- 
tively, given by: 



Iff) a = COSy| T)a + Sin y| |) Q 
\e) a = Silly | t) Q -COSy| i) a , 



(5) 



where the angle 9 a is given by the criterion tan# Q = 
A a /e a . We take the interaction Hamiltonian between 
the qubit and the TLS to be of the form: 



Hi 



X 



al q W TLS ) 



(6) 



III. INTUITIVE PICTURE 

We start our analysis of the problem by presenting a 
few physical principles that prove very helpful in intu- 
itively predicting the behaviour of the above-described 
system. Note that the arguments given in this section 
are well known [221 l23l l24j . For the sake of clarity, how- 
ever, we present them explicitly and discuss their roles 
in the problem at hand. 



A. Energy levels and eigenstates 

The first element that one needs to consider is the en- 
ergy levels of the combined qubit+TLS system. In order 
for a given experimental sample to function as a qubit, 
the qubit-TLS coupling strength A must be much smaller 
than the energy splitting of the qubit E q . We there- 
fore take that limit, as well as the limit A <C £tls> an d 
straightforwardly find the energy levels to be given by: 



£1 



Etls + E q A c 



E 2 = -y(ETL S -E q f + Xl + ^ 



1 



E 3 = +^\/(E TLS -E q r + Xl + ^ 



E x 



-Btls + E q A cc 
2 2 ' 



(7) 



where A cc = A cos q cos #tls , A ss = A sin q sin #tls ■ The 
corresponding eigenstates are given by: 



where A is the (uncontrollable) coupling strength between 
the qubit and the TLS. Note that any interaction Hamil- 
tonian that is a product of a qubit observable times a 



|i> = \gg) 

I 2 ) = cos||e 5 ) +sin^|#e) 
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|3) = sin-|ep) -cos-|#e) 

|4) = lee), " (8) 

where the first symbol refers to the qubit state and the 
second one refers to the TLS state in their respective 
uncoupled bases, the angle ip is given by the criterion 
tamp = A ss / (-Btls — E q ), and for definiteness in the 
form of the states |2) and |3) we have assumed that 
Etls > E q . 

Note that the mean-field shift of the qubit resonance 
frequency, A cc , is present regardless of the values of the 
qubit and TLS energy splittings. The avoided-crossing 
structure involving states |2) and |3), however, is only 
relevant when the qubit and TLS energies are almost 
equal. One can therefore use spectroscopy of the four- 
level structure to experimentally measure the TLS energy 
splitting -Etls and angle 6*tlSj as will be discussed in 
more detail in Sec. V. 



B. Rabi oscillations 

If a two-level system (e.g. a qubit) with energy split- 
ting ujq, initially in its ground state, is driven by a 
harmonically-oscillating weak field with a frequency to 
close to its energy splitting (up to a factor of K) as de- 
scribed by Eq. J3J), its probability to be found in the 
excited state at a later time t is given by: 



e ^ + (w _ Wo)2 2 W 

where = i/fig + (uj — ojq) 2 , and the on-resonance Rabi 
frequency Q = F\ sin(0 f - q )| (we take K = 1) 0. We 
therefore see that maximum oscillations with full g «-> e 
conversion probability are obtained when the driving is 
resonant with the qubit energy splitting. We also see 
that the width of the Rabi peak in the frequency domain 
is given by f2 - Simple Rabi oscillations can also be ob- 
served in a multi-level system if the driving frequency is 
on resonance with one of the relevant energy splittings 
but off resonance with all others. 



C. The effect of decoherence 

In an undriven system, the effect of decoherence is 
to push the density matrix describing the system to- 
wards its thermal-equilibrium value with time scales 
given by the characteristic dephasing and relaxation 
times. The effects of decoherence, especially dephasing, 
can be thought of in terms of a broadening of the energy 
levels. In particular, if the energy separation between the 
states |2) and |3) of Sec. III-A is smaller than the typical 
decoherence rates in the problem, any effect related to 
that energy separation becomes unobservable. Alterna- 
tively, one could say that only processes that occur on 



a time scale faster than the decoherence times can be 
observed. 

It is worth taking a moment to look in some more 
detail at the problem of a resonantly-driven qubit coupled 
to a dissipative environment, which is usually studied 
under the name of Bloch-equations [24l |25| . If the Rabi 
frequency is much smaller than the decoherence rates, the 
qubit will remain in its thermal equilibrium state, since 
any deviations from that state caused by the driving field 
will be dissipated immediately. If, on the other hand, the 
Rabi frequency is much larger than the decoherence rates, 
the system will perform damped Rabi oscillations, and it 
will end up close to the maximally mixed state in which 
both states \g) and |e) have equal occupation probability. 
In that case, one could say that decoherence succeeds in 
making us lose track of the quantum state of the qubit 
but fails to dissipate the energy of the qubit, since more 
energy will always be available from the driving field. 



D. Combined picture 

We now take the three elements presented above and 
combine them to obtain a simple intuitive picture of the 
problem at hand. 

Let us for a moment neglect the effects of decoherence 
and only consider the case w ~ ujq. The driving field 
tries to flip the state of the qubit alone. However, two of 
the relevant eigenstates are entangled states, namely |2) 
and 1 3). One can therefore expect that if the width of 
the Rabi peak, or in other words the on-resonance Rabi 
frequency, is much larger than the energy separation be- 
tween the states |2) and |3), the qubit will start oscillat- 
ing much faster than the TLS can respond, and the ini- 
tial dynamics will look similar to that of the uncoupled 
system. Only after many oscillations and a time of the 
order of (E3 — E2)" 1 will one start to sec the effects of the 
qubit-TLS interaction. If, on the other hand, the Rabi 
frequency is much smaller than the energy separation be- 
tween the states |2) and |3), the driving field can excite at 
most one of those two states, depending on the driving 
frequency. In that case the qubit-TLS interactions are 
strong enough that the TLS can follow adiabatically the 
time evolution of the qubit. In the intermediate region, 
one expects that if the driving frequency is closer to one 
of the two transition frequencies E2 — E\ and E3 — E\ , 
beating behaviour will be seen right from the beginning. 
If one looks at the Rabi peak in the frequency domain, 
e.g. by plotting the maximum g <-> e qubit-state con- 
version probability as a function of frequency, the single 
peak of the weak-coupling limit separates into two peaks 
as the qubit-TLS coupling strength becomes comparable 
to and exceeds the on-resonance Rabi frequency. 

We do not expect weak to moderate levels of decoher- 
ence to cause any qualitative changes in the qubit dy- 
namics other than, for example, imposing a decaying en- 
velope on the qubit excitation probability. As mentioned 
above, features that are narrower (in frequency) than the 
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decoherence rates will be suppressed the most. Note that 
if the TLS decoherence rates are large enough |2jj, the 
TLS can be neglected and one recovers the single Rabi 
peak with a height determined by the qubit decoherence 
rates alone. 

IV. NUMERICAL RESULTS 

In the absence of decoherence, we find it easiest to treat 
the problem at hand using the dressed-state picture • 
In that picture one thinks of the driving field mode as 



being quantized, and processes are described as involving 
the absorption and emission of quantized photons by the 
qubit+TLS combined system. As a representative case, 
which also happens to be the case of most interest to 
us, we take the frequency of the driving field to be close 
to the qubit and TLS energy splittings. For simplicity, 
we take those to be equal. We shall come back to the 
general case later in this section. Without going over the 
rather simple details of the derivation, we show the four 
relevant energy levels and the possible transitions in Fig. 
1. The effective Hamiltonian describing the dynamics 
within those four levels is given by: 



/ 



Q' -5uj + A cc - A ss /2 



V o 
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(10) 
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FIG. 1: Left: Energy levels and direct transitions between 
them in the dressed-state picture. Right: Energy splittings 
separated according to their physical origin. 

where Q' = Slo /2 3 / 2 , fio is the on-resonance Rabi fre- 
quency in the absence of qubit-TLS coupling, Su> — 
lj — ljq, and the Hamiltonian is expressed in the basis of 
states (|1, N), |2, N- 1), |3, N-l), |4, N-2}), where N is 
the number of photons in the driving field. We take the 
low temperature limit, which means that we can take the 
initial state to be |1, AT) without the need for any extra 
initialization. We can now evolve the system numerically 
and analyze the dynamics. After we find the density ma- 
trix of the combined qubit+TLS system as a function of 
time, we can look at the dynamics of the combined sys- 
tem or that of the two subsystems separately, depending 
on which one provides more insightful information. 

We start by demonstrating the separation of the Rabi 
peak into two peaks as the qubit-TLS coupling strength 
is increased. As a quantifier of the amplitude of Rabi os- 
cillations, we use the maximum probability for the qubit 
to be found in the excited state between times t = and 
t = 207r/17o, and we refer to that quantity as P^^ax- 
In Fig. 2 we plot -P|^ ax as a function of renormalized 
detuning Suj/Qo. As was explained in Sec. Ill, the peak 
separates into two when the qubit-TLS coupling strength 
exceeds the on-resonance Rabi frequency, up to simple 
factors of order one. The system also behaves according 
to the explanation given in Sec. Ill in the weak- and 



strong-coupling limits. When A is substantially smaller 
than f2oi oscillations in the qubit state occur on a time 
scale f2 , where f2 is the Rabi frequency defined in Sec. 
III-B, whereas the beating behaviour occurs on a time 
scale (£3 — -E^) -1 - When A is more than an order of mag- 
nitude smaller than fio, the effects of the TLS are hardly 
visible in the qubit dynamics within the time given above. 
On the other hand, when A is large enough such that the 
energy difference £3 — Ei is several times larger than fioj 
the dynamics corresponds to exciting at most one of the 
two eigenstates |2) and |3). We generally see that beating 
behaviour becomes less pronounced when the driving fre- 
quency is equal to the qubit energy splitting including the 



TLS-mean field shift, i.e. when lu = y Aq + (e q + A cc ) 2 , 

which corresponds to the top of the unsplit single peak 
or the midpoint between the two separated peaks. 

We also see some interesting features in the peak struc- 
ture of Fig. 2 that were not discussed in Sec. III. In the 
intermediate-coupling regime (Figs. 2b, c), we see a peak 
that reaches unit height, i.e. a peak that corresponds to 
full g <-> e conversion in the qubit dynamics at 5lo = 0. 
The asymmetry between the two main peaks in Fig. 2, 
as well as the additional dips in the double-peak struc- 
ture, were also not immediately obvious from the simple 
arguments of Sec. III. In order to give a first explanation 
of the above features, we plot in Fig. 3 a curve similar 
to that in Fig. 2(b) (with different 6>tls)> along with 
the same quantity plotted when the eigenstate |4, N — 2) 
is neglected, i.e. by using a reduced 3x3 Hamiltonian 
where the fourth row and column are removed from H e g . 
In the three-state calculation, there is no 5uj — peak, 
the two main peaks are symmetric, but we still see some 
dips. We also plot in Fig. 4 the qubit excitation probabil- 
ity as a function of time for the four frequencies marked 
by vertical dashed lines in Fig. 2(b). 



■5 




5co/n Q " 8co/n Q 

FIG. 2: Maximum qubit excitation probability P^ax be- 
tween t = and t = 207r/fi for A/fi = 0.5 (a), 2 (bj, 2.5 (c) 
and 5 (d). # q = 7r/4, and #tls = 7r/6. 



ingly narrow in that limit. In other words, the virtual 
intermediate state after the absorption of one photon is 
far enough in energy from the states 2} and |3) to make 
the peak invisibly narrow. It is rather surprising, how- 
ever, that in the intermediate-coupling regime the peak 
reaches unit g <-> e conversion probability, even though 
the transitions to states 2} and |3) are real, rather than 
being virtual transitions whose role is merely to mediate 
the coupling between states |1) and |4). We have ver- 
ified that the (almost) unit height of the peak is quite 
robust against changes in the angles 9 q and #tls f° r a 
wide range in A, even when that peak coincides with the 
top of one of the two main peaks. In fact, the Hamilto- 
nian H e g can be diagonalized rather straightforwardly in 
the case 5oj = 0, and one can see that there is no symme- 
try that requires full conversion between the states |1, N) 
and |4, N — 2). The lack of any special relations between 
the energy differences in the eigenvalues of H e s, however, 
suggests that almost full conversion should be achieved 
in a reasonable amount of time. 
























- 2 5co/n o 2 4 

FIG. 3: Maximum qubit excitation probability P^^ax be- 
tween t = and t = 20-k/Q.o for the four- level system 
(solid line) and the reduced three-level system (dashed line). 
A/fi<, = 2, 6 q = tt/4, and 6»tls = vr/5. 




FIG. 4: Qubit excitation probability as a function of 
time (solid line) for Slj/Q = (a), 0.92 (b), 1.61 (c) and 1.77 
(d). The dashed line is the occupation probability of state 
|4, N - 2>. A/n = 2, 9 q = tt/4, and 9 TLS = tt/5. 



By looking at Fig. 1, one might say that the Su> = 
peak clearly corresponds to a two-photon process cou- 
pling states |1) and 4). In fact, for further demonstra- 
tion that this is the case, we have included in Fig. 4 
the probability of the combined qubit+TLS system to 
be in state 4). This peak is easiest to observe in the in- 
termediate coupling regime. In the weak-coupling limit, 
the qubit and TLS are essentially decoupled, especially 
on the time scale of qubit dynamics. In the strong cou- 
pling limit, one can argue that a Raman transition will 
give rise to that peak. However, noting that the width 
of that peak is of the order of the smaller of the values 
£7q/A ss and £1q/A cc , we can see that it becomes increas- 



The asymmetry between the two main peaks in Figs. 
2 and 3 can also be explained by the fact that in one 
of those peaks state |4) is also involved in the dynamics 
and it increases the quantity -P-[ q ^ lax . As above, we have 
included in Figs. 4(b) and 4(c) the probability of the 
combined qubit+TLS system to be in state |4). 

In order to explain the dips in Figs. 2 and 3, we note 
that the plotted quantity, P^ aK , is the sum of four terms 
(in the reduced three-level system): a constant and three 
oscillating terms. The frequencies of those terms corre- 
spond to the energy differences in the diagonalized 3x3 
Hamiltonian. The dips occur at frequencies where the 
two largest frequencies are integer multiples of the small- 
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est one. Away from any such point, P^ax. w ^ u reacri 
a value equal to the sum of the amplitudes of the four 
terms. Exactly at those points, however, such a con- 
structive buildup of amplitudes is not always possible, 
and a dip is generally obtained. The width of that dip 
decreases and vanishes asymptotically as we increase the 
simulation time, although the depth remains unaffected. 

We also studied the case where the qubit and TLS en- 
ergy splittings were different. As can be expected, the 
effects of the TLS decrease as it moves away from reso- 
nance with the qubit. That is most clearly reflected in 
the two-peak structure, where one of the two main peaks 
becomes substantially smaller than the other. The two- 
photon peak was still clearly observable in plots corre- 
sponding to the same quantity plotted in Fig. 2, i.e. plots 
of -Pf^ax vs. d(j/Qo, even when the detuning between 
the qubit and the TLS was a few times larger than the 
coupling strength and the on-resonance Rabi frequency. 




-2-101234 



FIG. 5: Maximum qubit excitation probability P^^ax be- 
tween t = and t = 207r/f2o- The solid line corresponds to 
the case of no decoherence. The dotted (T^ = O.lfio /2n and 
r< T 2 LS) = 0.2n /27r), dashed (rj$ = and r^ T 2 LS) = 2fi ) 

and dash-dotted (F^ = S7o and r^ LS ' = 0) lines correspond 
to different decoherence regimes. A/Qq — 2, 6 q = 7r/4, and 
#tls = 7r/5. and are the relaxation and dephasing 

rates of the subsystem a, respectively. 

The truncated dressed-state picture with four energy 
levels is insufficient to study the effects of decoherence. 
For example, relaxation from state 2) to |1) does not 
necessarily have to involve emission of a photon into the 
driving-field mode. We therefore study the effects of de- 
coherence by treating the driving field classically. We 
then solve a Bloch-Redfield master equation with a time- 
dependent Hamiltonian and externally-imposed dephas- 
ing and relaxation times, as was done in Ref. [20(. In 
Fig. 5 we reproduce the four-level results of Fig. 3, i.e. 

max vs - S^/^o w ith no decoherence, along with the 
same quantity obtained when we take into account the ef- 



fects of decoherence. For a moderate level of decoherence, 
we see that the qubit excitation probability is somewhat 
reduced and all the features that are narrower than the 
decoherence rates are suppressed partially or completely 
by the effects of decoherence. For large qubit decoherence 
rates, the qubit excitation probability is greatly reduced 
close to resonance, where the Rabi frequency f2 takes its 
lowest values. The shallow dip in the dash-dotted line in 
Fig. (5) occurs because for those frequencies and in the 
absence of decoherence the maximum amplitude is only 
reached after several oscillations, whereas it is reached 
during the first few oscillations outside that region. For 
large TLS decoherence rates, the TLS becomes weakly 
coupled to the qubit, and a single peak is recovered in 
the qubit dynamics (with a height larger than either the 
two split peaks). All of these effects are in agreement 
with the simple picture presented in Sec. III. 

V. EXPERIMENTAL CONSIDERATIONS 

In the early experiments on phase qubits coupled to 
TLSs the qubit relaxation rate r[ q) (~40MHz) 

was comparable to the splitting between the two Rabi 
peaks A ss (~20-70MHz), whereas the on-resonance Rabi 
frequency f^o was tunable from 30MHz to 400MHz (note 
that, as discussed in Sec. Ill, the Rabi frequency can- 
not be reduced to values much lower than the decoher- 
ence rates, or Rabi oscillations would disappear alto- 
gether). The large relaxation rates in those experiments 
would make several effects discussed in this paper unob- 
servable. The constraint that f^o could not be reduced 
below 30MHz made the strong-coupling regime, where 
0,0 <C A ss , inaccessible. The weak-coupling regime, where 
fio 3> A ss , was easily accessible in those experiments. 
However, as can be seen from Fig. 2, it shows only a 
minor signature of the TLS. Although the intermediate- 
coupling regime was also accessible, as evidenced by the 
observation of the splitting of the Rabi peak into two 
peaks, observation of the two-photon process and the 
additional dips of Fig. 2 discussed above would have re- 
quired a time at least comparable to the qubit relaxation 
time. That would have made them difficult to distinguish 
from experimental fluctuations. 

With the new qubit design of Ref. 26], the qubit re- 
laxation time has been increased by a factor of 20. The 
constraint that Qq must be at least comparable to T± % 
no longer prevents accessibility of the strong-coupling 
regime. Furthermore, since our simulations were run for a 
period of time corresponding to approximately ten Rabi 
oscillation cycles, i.e. shorter than the relaxation time 
observed in that experiment, all the effects that were 
discussed above should be observable, including the ob- 
servation of the two-photon peak and the transition from 
the weak- to the strong-coupling regimes by varying the 
driving amplitude. 

We finally consider one possible application of our re- 
sults to experiments on phase qubits, namely the problem 
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of characterizing the environment composed of TLSs. As 
we shall show shortly, characterizing the TLS parame- 
ters and the nature of the qubit-TLS coupling are not 
independent questions. The energy splitting of a given 
TLS, -EtlSj can be obtained easily from the location of 
the qubit-TLS resonance as the qubit energy splitting is 
varied. One can then obtain the distribution of values of 
-Etls for a large number of TLSs, as was in fact done in 
Ref. |2(j . The splitting of the Rabi resonance peak into 
two peaks by itself, however, is insufficient to determine 
the values of Atls an d £tls separately. By observing 
the location of the two-photon peak, in addition to the 
locations of the two main peaks, one would be able to de- 
termine both A cc and A ss for a given TLS, as can be seen 
from Fig. 1. Those values can then be used to calculate 
both -Etls and #tls = arctan(ATLs/ £ TLs) of that TLS. 
The distribution of values of #tls can then be used to 
test models of the environment, such as the one given in 
Ref. U3 to describe the results of Ref. [27|. 

In order to reach the above conclusion, we have made 
the assumption that the distribution of values of $tls for 
those TLSs with sufficiently strong coupling to the qubit 
is representative of all TLSs. Since it is generally be- 
lieved that strong coupling is a result of proximity to the 
junction, the above assumption is quite plausible, as long 
as the other TLSs share the same nature. Although it is 
possible that there might be two different types of TLSs 
of different nature in a qubit's environment, identifying 
that possibility would also be helpful in understanding 
the nature of the environment. We have also assumed 
that 9 q does not take the special value ir/2 (note that, 
based on the arguments of Refs. 0, |23, we are also 
assuming that generally #tls t/2). That assumption 
would not raise any concern when dealing with charge 
or flux qubits, where both A q and e q can be adjusted 
in a single experiment, provided an appropriate design 
is used. However, the situation is trickier with phase 
qubits. The results in that case depend on the nature of 
the qubit-TLS coupling, which we discuss next. 

The two mechanisms that are currently considered the 
most likely candidates to describe the qubit-TLS cou- 
pling are through either (1) a dependence of the Joseph- 
son junction's critical current on the TLS state or (2) 
Coulomb interactions between a charged TLS and the 
charge across the junction. In the former case, one has 
an effective value of 9 q that is different from n/2 [2^ , and 
the assumption of an intermediate value of 9 q is justified. 
In the case of coupling through Coulomb interactions, on 
the other hand, one effectively has 6 q — tt/2, and there- 
fore A cc vanishes for all the TLSs. In that case the two- 



photon peak would always appear at the midpoint (to a 
good approximation) between the two main Rabi peaks. 
Although that would prevent the determination of the 
distribution of values of #tls j it would be a strong indi- 
cation that Coulomb interactions with the charge across 
the junction are responsible for the qubit-TLS coupling 
rather than the critical current dependence on the TLS 
state. Note also that if it turns out that this is in fact 
the case, and the distribution of values of #tls cannot 
be extracted from the experimental results, that distri- 
bution might be irrelevant to the question of decoherence 
in phase qubits. 



VI. CONCLUSION 

We have studied the problem of a harmonically-driven 
qubit that is interacting with an uncontrollable two-level 
system and a background environment. We have pre- 
sented a simple picture to understand the majority of 
the phenomena that are observed in this system. That 
picture is composed of three elements: (1) the four-level 
energy spectrum of the qubit+TLS system, (2) the ba- 
sic properties of the Rabi-oscillation dynamics and (3) 
the basic effects of decoherence. We have confirmed the 
predictions of that picture using a systematic numeri- 
cal analysis where we have varied a number of relevant 
parameters. We have also found unexpected features 
in the resonance-peak structure. We have analyzed the 
behaviour of the system and provided simple explana- 
tions in those cases as well. Our results can be tested 
with available experimental systems. Furthermore, they 
can be used in experimental attempts to characterize the 
TLSs surrounding a qubit, which can then be used as part 
of possible techniques to eliminate the TLSs' detrimental 
effects on the qubit operation. 
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